Abstract. Results of a 2000 August 12-13 BeppoSAX 
observation of the 7.1 hr eclipsing, dipping, bursting, tran- 
sient, low-mass X-ray binary (LMXRB) X 1658-298 are 
presented. The spectrum outside of eclipses, dips and 
bursts can be modeled by the combination of a soft disk- 
blackbody and a harder Comptonized component with a 
small amount (1.3 x 10^^ atom cm~^) of low-energy ab- 
sorption. In contrast, an RXTE observation 18 months 
earlier during the same outburst, measured an absorp- 
tion of 5.0 X 10^^ atom cm~^. Such a change is consis- 
,—1 . tent with a thinning of the accretion disk as the out- 
^^ burst progresses. Structured residuals from the best-fit 
^^ i; spectral model are present which are tentatively identi- 
^^ ' fied with Ne-K/Fe-L and Fe-K shell emission. The spec- 
'^ , tral changes during dips are complex and may be modeled 
1-^ ' by a strong (~3 x 10^'^ atom cm~^) increase in absorption 
QQ ■ of the Comptonized component only, together with reduc- 
,— I tions in normalizations of both spectral components. This 
■ behavior is in contrast to the "complex continuum" model 
for X-ray dip sources, where the softer blackbody com- 
ponent rapidly suffers strong absorption. It is however, 
X^ similar to that found during recent XMM-Newton obser- 
Cn ' vations of the eclipsing, dipping, LMXRB EXO 0748-676. 
t^^ ' An updated orbital ephemeris is provided which does not 
^^ , require a quadratic term, in contrast to that of Wachter 
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1. Introduction 

X 1658-298 is a 
SAS-3 (Lewin et al 

X 1658-298 in 1976 October and 1977 June 
was a pparently no significant persistent emission (Lewin 
During 1978 March, observations with SAS-3 and 



transient X-ray source discovered by 
19761 ). Bursts were reported from 

but there 



197C) 



HEAO-1 detected a clearly modulated persistent emission 
which showed erratic intensity variations followed by 15 
minute duration X-ray eclipses which repeated periodi- 
cally (Cominsky & Wood 1984). By combining data span- 
ning ~2 years an orbital period of 7.11 h was obtained 
(Cominsky & Wood 1989|) . Following the discovery using 
EXOSAT of a second bursting, di pping , eclipsing, tran- 
sient (EXO 0748-676; Parmar et al. |1986D the region of sky 
containing X 1658-298 was observed by EXOSAT. How- 
ever, the source was not detected. Indeed, until recently 
X 1658-298 has not been observed to be active. In 1991 
a 25 ks ROSAT Position Sensitive Proportional Counter 
observation failed to detect X 1658-298 with a 0.2-2.4 keV 
upper limit of ^2 x 10"^'* erg cm~^ s~^. However, on 1999 
April 2, after a 21 year quiescent interval, an 18 mCrab 
(2-9 keV) source at a position consistent with X 1658-298 
was detected in one of the BeppoSAX Wide Field Cameras 
(in 't Zand et al. 1999a). Observations a day later con- 



firmed that the V=18.3 optical counterpart to X 1658-298 



(V2134 Oph; Doxsey et al. |1979D had brightened signifi- 
cantly and was exhibiting high excitation He ii emission, 
confirming that t he so urce was undergoing an outburst 
(Augusteijn et al. |l999|) . 

X 1658-298 was observed by RXTE a few days later 
(Wachter et al. 2000, hereafter WSB) confirming that the 
source had resumed its strong and persistent X-ray emis- 
sion. WSB combined their RXTE eclipse timing measure- 
ments with those obtained with HEAO A-1 and SAS-3 and 
found that the orbital period of X 1658-298 had decreased 
with a timescale of 10^ yr. The eclipse transitions of 6 
to 13 s were clearly resolved in the RXTE measurements. 
Contemporaneous observations revealed that the optical 
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modulation is highly variable, and exhibits a 0.2 mag- 
nitude eclipse feature coincident with the X-ray eclipse. 
Recently, nearly coherent pulsations were discovered from 
this s ource during Type I X-ray bursts (Wijnands et al. 
^000 ). These imply a spin period of ~ 1.8 ms, if the oscil- 
lations are directly related to the neutron star rotation. 

Only two Low-Mass X-ray Binary (LMXRB) systems 
are known that exhibit both dips and eclipses. (About 
10 LMXRBs exhibit only dips). During dips the observed 
spectral changes are complex and inconsistent with simple 
absorption from cold material, as might be expected. De- 
tailed modeling of these spectral changes provides a pow- 
erful means of studying the structure and location of the 
emitting and absorbing regions in LMXRB (e.g., Parmar 



et al. 1986; Church & Balucihska-Church 1995; Church et 
al. 1997). This modeling has revealed that at least two 



emission components are required. These regions appear 
to undergo different amounts of absorption during dips 
with a point-like blackbody undergoing rapid obscura- 
tion and a cutoff power-law component undergoing a more 
gradual increase in obscuration during dipping intervals. 
This gradual change is highly suggestive of an extended 
component. We refer to this as the complex-continuum 
model. The same spectral model has been used to describe 
the ASCA spectra of a large number of LMXRBs (Church 
& Balucinska-Church p00l| ). 

In the case of EXO 0748-676, (the other eclipsing and 
dipping LMXRB), recent XMM-Newton Reflection Grat- 
ing Spectrometer (RGS) observations have revealed the 
presence of narrow spectral features due to ionized neon, 



oxygen and nitrogen (Cottam et al. 2001). Simultaneous 
European Photon Imaging Camera (EPIC) data revealed 
that the 0.2-10 keV spectra can be fit with a two compo- 
nent model consisting of a central Comptonizcd compo- 
nent and a more extended thermal halo (Bonnet-Bidaud 



et al. 2001). In contrast to the complex continuum contin- 



uum model, the spectral variations observed during dips 
can be mainly accounted for by variations in the absorbing 
column affecting the Comptonized component only. 
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Fig. 1. 2-10 keV MECS lightcurve of X 1658-298 with a binning of 32 s. 2 completely and 2 partially observed eclipses, 
irregular intensity dips and 9 X-ray bursts are evident. The peak heights of the bursts are truncated. The arrows mark 
the expected mid-eclipse times. Time is hours of 2000 August 12 
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Fig. 2. The folded hardness ratio (top panel) and 2-10 keV MECS lightcurve showing the deep dipping activity 
between phases ~0. 68-0. 95 and the narrow eclipse centered on phase 0. X-ray bursts have bee n rem oved prior to the 
folding. The horizontal lines indicate phases used to extract the three dip-spectra (see Sect. 5.2). The central line 
corresponds to the deep dip interval, the two middle lines to the medium dipping, and the two outer lines to the 
shallow dipping interval 



2. Observations 

Results from the Low-Energy Concentrator Spectrometer 
(LECS; 0.1-10 keV; Parmar et al. |1997D , the Medium- 
Energy Concentrator Spectrometer (MECS; 1.8-10 keV; 
Boella et al. 1997), the High Pressure Gas Scintillation 
Pr oport ional Counter (HPGSPC; 5-120 keV; Manzo et 
al. |1997| ) and the Phoswi ch De tection System (PDS; 15- 



300 keV; Frontera et al. |1997| ) on-board BeppoSAX are 



presented. All these instruments are coaligned and col- 
lectively referred to as the Narrow Field Instruments, or 
NFL The MECS consists of two grazing incidence tele- 
scopes with imaging gas scintillation proportional coun- 
ters in their focal planes. The LECS uses an identical con- 
centrator system as the MECS, but utilizes an ultra-thin 
entrance window to extend the response down to 0.1 keV. 
The non-imaging HPGSPC consists of a single unit with 
a collimator that was pointed on-source during the whole 
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observation. The non-imaging PDS consists of four inde- 
pendent units arranged in pairs each having a separate 
coUimator. Each colhmator was alternatively rocked on- 
and 210' off-source every 96 s during the observation. 

The region of sky containing X 1658-298 was observed 
by BeppoSAX on 2000 August 12 10:33 UT to August 13 
18:14 UT. Good data were selected in the standard way 
using the SAXDAS 2.0.0 data analysis package. LEGS and 
MEGS data were extracted centered on the (on-axis) po- 
sition of X 1658-298 using radii of 8' and 4', respectively. 
The exposure times in the LEGS, MEGS, HPGSPG, and 
PDS instruments are 20.2 ks, 53.5 ks, 50.1 ks, and 24.3 ks, 
respectively. Background subtraction for the imaging in- 
struments was performed using standard files, but is not 
critical for such a bright source. Background subtraction 
for the HPGSPG was carried out using data obtained 
when the instrument was observing the dark Earth and 
for the PDS using data obtained when the collimator was 
offset from the source. 
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Fig. 3. The X 1658-298 eclipse timing residuals with re- 
spect to the best-fit linear ephemeris. The BeppoSAX 
measurements are the two most recent data points. The 
best-fit quadratic ephemeris (long dashed line) is also in- 
dicated 



3. X-ray lightcurve 

Fig. I shows the 2-10 keV MEGS hghtcurve of X 1658-298 
with a binning of 32 s spanning just over 4 orbital cycles. 
As expected, eclipses, irregular dips in X-ray intensity, 
and bursts are present. Neither the dips nor the eclipses 
are quite total. A folded 2-10 keV lightcurve with bursts 
removed, and the (4.0-10.0 keV/2. 0-4.0 kcV) hardness ra- 
tio as a function of phase are plotted in Fig. 0. From this 
figure it is evident that strong dipping behavior occurs 
between phases 0.68-0.95 (where center of X-ray eclipse 
corresponds to phase 0.0). At other phases, the hardness 
ratio increases slowly with increasing phase, from a mini- 
mum around phase 0.0. This may indicate that low-level 
dipping is present throughout the orbital cycle. Alterna- 
tively, this variation could be due to a dependence of the 
X-ray emission on orbital phase. 

4. Eclipse timing 

The arrival times of the eclipses were determined by fitting 
a "step and ramp" function (see Parmar et al. 1986) to the 
MEGS light curves (which were obtained at 8 s resolution) 
of the two complete eclipses. This model has 7 free param- 
eters: the start- and end-times of the ingress and egress, 
and the count rates before, during, and after eclipse. The 
time of eclipse was calculated as the mean between the 
end time of ingress and the start time of egress. 

The derived eclipse times, which extend the measure- 
ments of WSB by more than a year, or 1600 cycles, 
are summarized in Table nl These times were then com- 
bined with the eclipse times compiled by WSB and lin- 
ear and quadratic ephemerides investigated. Fig. H shows 
the residuals with respect to both ephemerides. Since the 
WSB results require a quadratic solution (see below) a 
quadratic ephemeris was first fit to the eclipse times to 



obtain a ^ of 45.4 for 5 degrees of freedom (dof). The up- 
dated quadratic ephemeris for the eclipse centers is given 
by: 
Toci = 2443059.22609(11) + 0.29650453(33) N 

+ (1.1 ±0.3) X lO^^^N^. 
Note that the sign of the quadratic term is opposite from 
that of WSB, and the magnitude a factor '^10 smaller. A 
linear ephemeris was next fit to the eclipse times to give 
a x^ of 62.9 for 6 dof: 

Tod = 2443059.225826(93) + 0.2965045746(34) N 
The x^ values of the two fits indicate that the difference in 
X^ is not significant (the probability of finding a decrease 
in the x^ by this value (or larger) is 22%). Note that the 
X^ values are much larger than 1, indicating that neither 
fit is formally acceptable. 

To investigate the large difference in the quadratic 
terms of the two ephemerides, only the data points re- 
ported in WSB were fit. This gives the same solution as 
WSB found (Fig. |). The WSB solution goes through the 
two sets of data points, but deviates strongly from the lin- 
ear ephemeris in the interval where no measurements are 
available. From this we conclude that the ephemeris ob- 
tained by WSB is an artifact of the strongly non-uniform 
sampling of the data and that their P has no physical rele- 
vance. In Fig. the linear ephemeris when the BeppoSAX 
data are included is plotted. From this we conclude that a 
linear ephemeris is an acceptable solution, although devi- 
ations (of the order of a few tr) are present. These may be 
similar to the "random" deviati ons in eclipse times seen 
from EXO 0748-676 (Hertz et al. |l997|) . 



5. X-ray spectrum 
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Fig. 4. The X 1658-298 eclipse timing residuals with re- 
spect to a linear ephemeris obtained using the WSB data 
only. The WSB quadratic ephemeris is indicated by the 
long-dashed line. The BeppoSAX data points are strongly 
inconsistent with the WSB quadratic ephemeris. Note the 
different vertical scale with respect to Fig. y 

Table 1. New X-ray center of eclipse times (corrected 
to the solar system barycenter) for X 1658-298. A compi- 
lation of previous measurements can be found in WSB. 
Cycle numbers are given using the convention in WSB. 
Uncertainties are given at 68% confidence 



JDfT 



Uncertainty Cycle 



2,451,769.344257 
2,451,769.937259 



0.000016 
0.000015 



29376 
29378 



5.1. Overall spectrum 

The overall spectrum of X 1658-298 was first investigated 
by simultaneously fitting data from all the BeppoSAX 
NFL Only intervals when the source was not bursting, dip- 
ping, or eclipsing were included. This resulted in exposure 
times in the LEGS, MEGS, HPGSPG, and PDS instru- 
ments of 13.0, 30.8, 27.6, and 13.5 ks, respectively. In or- 
der to ensure applicability of the x^ statistic, the extracted 
spectra were rebinned such that at least 20 counts per bin 
were present and such that the resolution of the instru- 
ment is not oversampled by more than a factor 3. LEGS 
data were used in the energy range 0.5-4.0 keV, MEGS 
data between 1.85-10.5 keV, HPGSPG data between 8-20 
keV, and PDS data in the 15-30 keV energy range, where 
the instrument responses are well determined and suffi- 
cient counts obtained. This gives background-subtracted 
count rates of 3.5, 6.8, 3.1 and 0.7 s~^ in these energy 
ranges for the LEGS, MEGS, HPGSPG, and PDS, respec- 
tively. The photo-electric absorption cross sections of Mor- 



rison & McGammon ( 1983| ) were used throughout. Factors 
were included in the spectral fitting to allow for normaliza- 
tion uncertainties between the instruments. These factors 
were constrained to be within their usual ranges during 
the fitting. All spectral uncertainties and upper-limits are 
given at 90% confidence. A systematic uncertainty of 1% 
was added quadratically to account for known instrument 
calibration uncertainties. 

Since the source is 7° away from the galactic plane 
we verified that e mission from the galactic ridge (e.g., 
Valinia & Marshall 199S) is unlikely to affect the back- 



ground determination in the non-imaging instruments. 
The HPGSPG background data are obtained during ob- 
servations of the dark Earth and therefore the resultant 
spectrum will be the sum of the source spectrum and the 
galactic ridge spectrum. However, given the large offset 
from the galactic plane and the brightness of the source, 
this is unlikely to be a significant effect. The PDS back- 
ground spectra are obtained by rocking the instrument 
collimator 3?5 in opposite directions. Thus, changes in the 
galactic ridge emission (as a function of galactic latitude) 
could lead to an imperfect background subtraction. We 
therefore accumulated PDS spectra using only the posi- 
tive and negative offset background directions. The differ- 
ences in the PDS spectra obtained in this way are much 
smaller than the statistical uncertainties. 

The resulting spectra were first fit with simple 
models including absorbed power-laws with photon- 
index a, thermal bremsstrahlung and cutoff power-laws 
(E~"exp— (Ec/kT)). This showed that no such model 
provides an acceptable fit with e.g., an absorbed power- 
law with a ~ 2 giving a x^ of 1149 for 103 dof. A 
cutoff power-law with a = 1.3 and absorption, Njj, of 
3.5 X 10^^ atom cm"'^ provides a significantly better fit 
with a x^ of 231 for 102 dof implying the presence of a 
high-energy break at 6.1 keV. More complex models, con- 
sisting of combinations of soft and hard components were 
next tried. The soft components tried were blackbodies, 
coUisionally ionized thermal plasmas (the mekal model 
in XSPEC ), and the disk-blackbody model of Mitsuda et 
al. ( 1984 ). This model assumes that the gravitational en- 
ergy released by the accreting material is locally dissipated 
into blackbody radiation, that the accretion fiow is contin- 
uous throughout the disk, and that the effects of electron 
scattering are negligible. There are only two parameters 
of the model, rin(cosi)*''^ where r-m is the innermost ra- 
dius of the disk, i is the inclination angle of the disk and 
fcTin the blackbody effective temperature at Tin. The hard 
components tried were power-laws, cutoff power-laws and 
the XSPEC COMPTT Gomptonization model described in 
Titarchuk (1994). This model contains as free parame- 



ters the temperature of the Gomptonizing electrons fcTo, 
the plasma optical depth with respect to electron scat- 
tering To and the input temperature of the soft photon 
(Wien) distribution kTyj. A spherical geometry was as- 
sumed for the comptonizing region. Of all the possible 
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Fig. 5. The 0.5-30 keV non-dip, non-burst, and non-eclipse ("persistent emission") X 1658-298 count spectrum and 
residuals when a disk-blackbody and COMPTT model are fit. The right hand panel shows the deconvolved photon 
spectrum with the contributions of the two components indicated 



combinations, the best-fit two component model is a disk- 
blackbody and a COMPTT component with a x^ of 138.6 
for 99 dof (Fig. ^). The absorption corrected 2-10 keV 
luminosity is 1.6 x 10^^ erg s~^ for a distance of 15 kpc 



(Cominsky 1981[ ). Table H gives the fit results for mod- 
els with x^<2- The absorption derived using the best- 
fit model of (1.3 ± 0.2) x 10^^ atom cm~^ is similar to 
the galactic absorption in the direction of X 165 8-298 of 
1.8 X 10^1 atom cm'^ (Dickey & Lockman |l990| ) . 

We note here that we only use the COMPTT model 
as a convenient fit model. The physical interpretation of 
the obtained parameters (which we do not attempt to 
do) is difficult. The COMPTT model assumes the Wien- 
approximation for the input photons; for large values of 
kTo and low energies this breaks down. In our case the 
COMPTT model underestimates the flux up to 40% at the 



lowest energies (see Church & Balucihska-Church 2001). 
Additionally, there is evidence that the size of the Comp- 
toniz cd region is quite large (typically 50000 km; Church 
2001). The temperature of the seed-photons if they origi- 
nate from the disk at these distances is much lower than 
best-fit seed-temperature for typical luminosities. Addi- 
tionally, the difference in shape (at low energies) as dis- 
cussed above will influence the normalization and param- 
eters of the second (disk-blackbody) component, making 
a physical interpretation of these parameters even harder. 

Interestingly, recent results from the XMM-Newton 
EPIC show that the soft component seen from EXO 0748- 



676 may be modeled by a coUisionally ionized ther- 
mal plasma with strongly non-solar abundances (Bonnet- 
Bidaud et al. ^001] ). This is not the case in X 1658-298, 
where replacing the disk-blackbody by a ^^1.4 keV mekal 
component with abundances fixed at solar gives a signif- 
icantly worse fit with a x^ of 241.8 for 99 dof. Allowing 
the abundances to vary does not provide a signifi cantly 
improved fit. Additionally, Bonnet-Bidaud et al. ( 2001 ) 
find that their EXO 0748-676 power-law component suf- 
fers extra absorption of '^6 x 10^^ atom cm~^ during non- 
dipping or flaring intervals. No such additional absorption, 
with a 90% confidence upper limit of 9 x 10^° atom cm~^, 
is required by the X 1658-298 COMPTT component during 
non-dipping intervals. 

We next investigated the reasons for the poor (x^ = 
138.6 for 99 dof) best-fit. Inspection of Fig. || reveals 
the presence of structured residuals, with amplitudes of 
^^5% with similar forms for the different instruments in 
the overlapping energy ranges. This indicates that the ob- 
served inadequacies of the spectral models may be phys- 
ical, and not caused by incorrectly modeled instrumental 
responses. In order to investigate the nature of the resid- 
uals we selected spectra (i) based on count rate, and (ii) 
based on orbital phase. The first check was performed in 
order to see if the residuals are somehow caused by accu- 
mulating data from a range of source intensities. However, 
no difference in the structure of the residuals was found 
between the high and low count rate selected spectra. The 
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second test was performed since it appears that the hard- 
ness ratio increases close to the deep dipping phases (see 
Fig. H). Again, no difference in the structure of the resid- 
uals was found between spectra obtained at "early" and 
"late" phases. We therefore conclude that the observed 
residual structure is probably an intrinsic feature of the 
X 1658-298 spectrum. 

Finally, we tried the ad-hoc addition of Gaussian emis- 
sion features. Including narrow features at 1.06 ±0.08 and 
6.23 ± 0.13 keV with equivalent widths of 26 ± 16 and 
41 ± g'^ eV reduced the x^ to 114.7 for 95 dof. The proba- 
bility of obtaining a decrease in x^ by this value (or larger) 
is 1.2%. These features may be tentatively identified with 
Ne-K/Fe-L and Fe-K shell emission. Inspection of the re- 
maining residuals suggests that other discrete spectral fea- 
tures may be present in the spectrum. Given that XMM- 
Newton RGS observations have revealed O, Ne, and N 
features in the spectrum of the other eclipsing, dipping, 



'deep" dip spectrum 



LMXRB EXO 0748-676 (Cottam et al. |200l| ), this is per- 
haps unsurprising. We verified that adding similar Gaus- 
sians features to the other 3 models listed in Table does 
not produce markedly better fits. In all cases the change 
in X 
that of the best-fit model 



^ when the Gaussians were included was similar to 



5.2. Dip spectrum 

In order to understand the spectrum of a dip source, a suit- 
able model must give acceptable fits to both non-dip and 
dipping intervals. Since dips are almost certainly due to 
obscuration by intervening material, such a model should 
fit the dip spectra without requiring any of the parame- 
ters that characterize the shape of the emission, such as 
the temperature or power-law slope to be changed. These 
parameters are usually obtained by fitting to the non-dip 
spectrum. In the case of the "complex-continuum" spec- 
tral models, which have been used to model the spec- 
tral changes seen in many of the dip sources, the emis- 
sion consists of a point-like blackbody plus an extended 
Comptonized component. During dips the blackbody is 
rapidly absorbed, consistent with it being point-like, while 
the Comptonized component is more progressively ab- 
sorbed, consistent with i t being extend ed (e. g.. Church & 
Balucihska-Church |l99|; Church et al. [19971) . The black- 
body is usually assumed to originate from the surface of 
the neutron star or from a boundary layer between the 
neutron star and an accretion disk, and the extended com- 
ponent from an isothermal accretion disk corona. 

In order to investigate the spectral changes during 
X 1658-298 dipping intervals, three sets of LEGS and 
MEGS spectra were generated corresponding to shallow, 
medium, and deep dipping intervals (see Fig. O for a spec- 
trum during the deep dip). Selection was made using the 
orbital phases indicated in Fig. |. The MEGS 2-10 keV 
count rates were 5.4, 3.4, and 0.8 s~^ with exposures of 
5.2, 5.4, and 4.2 ks, respectively. Care was taken to min- 
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Fig. 6. Top panel: The LEGS and MEGS spectra of 
X 1658-298 obtained during the deep dip (see text). The fit 
consists of the model described in the text, without the ad- 
dition of the Fe-K line to clearly show the residuals near 
6.4 keV. Lower panel: the photon spectrum during the 
deep dip with the separate components indicated. Com- 
pare this with Fig. || to see the large effect of absorption 
on the COMPTT component. The Fe line is included in the 
model for this panel 



imize the range of count rates included in each accumu- 
lation, since this could lead to "mixing" of intervals with 
different spectra. Due to the low count rates and exposure 
times HPGSPC and PDS data were not included. 



5.2.1. Simple absorption 

Initially, the best-fit persistent emission model, modified 
to allow additional absorption to both components, was 
fit to the dip spectra. It was not possible to obtain satis- 
factory fits. Inspection of the residuals revealed that sig- 
nificant soft excesses were present in all 3 dip spectra. 
These can best be modeled if the absorptions (N^"'""'" and 
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Table 2. X 1658-298 fit results to the non-dip, non-eclipse, non-burst spectrum. Nh is in units of 10^^ atom cm~^. 
r'in(cosi)°'^ is in units of km for a distance of 15 kpc. 90% confidence limits are given. CPL = cutoff power-law, CTT= 
COMPTT, BB = blackbody, DBB = disk-blackbody. All temperatures are in keV. The normahzation ("XSPEC units") 
of the COMPTT component in the best-fit (DBB+CTT) model is (3.51 ± '^'l^) x 10"^ 



Model 



N 



H 



kTBB/kTln nn(c0Si)° 



Ec (keV) kT^ 



kT, 



xVdof 



BB-HCPL 
DBB+CPL 
BB+CTT 
DBB-HCTT 



0.29 ±0.03 
0.27 ±0.03 
0.07 ±0.02 
0.13 ±0.02 



0.77 ±0.07 
1.17 ± 
0.49 ± 
1.63 ±0.13 



0.15 
0.07 
0.01 
0.04 



3.7 ± 0.8 



2.3 ± 



0.5 
0.4 



1.25 ±0.07 
1.09 ±0.16 



7.1 ±0.9 

7.2 ± 



1.6 
1.0 



0.92±g:?f 
0.39 ±0.03 



4.8 ±0.9 

4.0 ± g:i 



167.1/100 

158.7/100 

3.5 ± 0.6 144.5/99 

5.4 ±0.6 138.6/99 



N§^^) and normalizations (/ctt and /dbb) of the disk- 
blackbody and COMPTT components are allowed to vary 
separately. This is the simplest model that is able to ad- 
equately model the spectral changes during dipping in- 
tervals and is referred to here as the "simple" absorption 
model. 

The results of fitting this simple absorption model are 
reported in Table ^ and show that during dips the nor- 
malizations of both components decrease, while only the 
COMPTT component is affected by strong additional ab- 
sorption. The depth of dipping is determined primarily by 
the reductions in normalization and not by an increase 
in absorption. In Fig. m the normalizations of both com- 
ponents are plotted as a function of the total flux in the 
dip spectra. This shows that the disk-blackbody normal- 
ization has a stronger dependence on total flux than the 
COMPTT component decreasing to a value of ~5% of the 
persistent emission during the deep dip. In contrast, the 
COMPTT component has a value ~60% of the persistent 
normalization during the deep dip. 

We investigated the sudden "onset" of the absorption 
(i.e. even the shallow-dip spectrum shows COMPTT ab- 
sorption of 3.4 X lO^'^ atom cm^^) by fitting the quiescent 
spectrum with the same simple absorption model. We were 
unable to obtain a fit to this spectrum with strong extra 
absorption for either component, but instead the best-fit 
was identical to that reported in Table |[ We also exper- 
imented with several other selections to obtain the dip 
spectra (using only 2 phase intervals, or using 4 intervals 
within the dip based on count rate). All these selections 
gave similar results to those described above. 

We also tried fitting the dip spectra with the same 
model, but with either the temperature of the COMPTT 
component (kTe) or the disk-blackbody temperature 
(A;Tin) as free parameters. This gives almost identical re- 
sults to those given in Table 0, with the possible exception 
of the fit to the deep-dip spectrum when kT^ is treated 
as a free parameter. This may be applicable if e.g., dif- 
ferent parts of a structured emitting region are covered 
during dipping. This gives a x^ of 61.7 for 51 dof with 
kTc — 4.6 keV, while the normahzation of the COMPTT 
component increases to 3.53x10"^, and N^"'""'" increases 
to 3.6 X 10^^ atom cm~^. However, the fact that in only 




2.0 3.0 

2-10 keV Flux (10 



Fig. 7. The normalizations of both the disk-blackbody 
(r?jCosi, left vertical scale, filled squares) and COMPTT 
components (right vertical scale, filled circles). From left 
to right the points correspond to the deep dip, medium 
dip, shallow dip, and persistent spectra. The rightmost 
circle has been shifted by a small amount for clarity. See 
text for a discussion on the value and uncertainty for the 
rightmost (i.e. persistent) disk-blackbody normalization. 
The points are connected with dotted lines to guide the 
eye 



1 out of 6 trials (three spectra with 2 models) the x^ is 
reduced with a significance of 6% probably means that 
there is no compelling reason to treat the temperature of 
either the COMPTT or disk-blackbody components as free 
parameters during dips. 

The uncertainty in the value of rin(cosi)*'-^ obtained 
from the fits to the dip spectra (see Table ph is much 
smaller than that obtained from the fits to the persis- 
tent emission spectrum (Table |^), although there are far 
fewer counts in the dip spectra. This results from fixing 
the two temperatures in the model which is used to fit the 
dip spectra, while the temperatures are free parameters in 
the fit to the persistent emission spectrum. The high value 
for Tin cos i obtained from the shallow dip spectrum (com- 
pared to that obtained from the quiescent spectrum), is 
hard to understand physically. However, there is a strong 
correlation between Tin and rin cosi and fixing Tin to 1.76 
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Table 3. Results of fitting the simple absorption model to the dip spectra, rin(cosi)'''^ assumes a distance of 15 kpc 



Parameter 



Deep 



Medium 



Shallow 



j^CTT (jo22 ^^^^ ^j^-2) 




30.8 ±4.7 


23.1±il 


33.6±^:i 


/CTT (10-2) 




2.i2±g:l! 


3.15 ±0.23 


3.49±g:i 


Ng™ (10^2 atom cm'^) 




<0.07 


<0.16 


0.07±0.04 


rin(cosi)«-5 




0.77±0.03 


9 91 1 0.03 


3.02±0.03 


Flux (2-10 keV) (erg cm"^ s" 


-') 


0.95 X 10"^° 


3.3 X 10"^° 


5.0 X 10"^° 


xVdof 




66.1/52 


107.2/69 


74.8/77 



Table 4. Fit results of the complex absorption model to 
the dip spectra. Units are the same as in Table ^ 



Parameter Deep 



Medium 



Shallow 



f 
xVdof 



37.3 ±2.5 
123 ± II 
0.91 ±0.01 
70.6/53 



17.8 ±2.3 
>148 

0.31 ±0.02 
115.5/70 



<0.09 

13.1 ± i:| 

0.47 ±0.04 
75.7/79 



(best-fit value ±90% confidence uncertainty) results in a 
value for rjn cos i which is within the uncertainty quoted 
in Table |. 

We note that the x^ for the medium dip spectrum is 
formally unacceptable with a 0.2% probability that the 
model describes the data. This might be caused by the 
large amount of spectral variability during the intervals 
over which the medium dip spectrum was selected (see 
Fig. d). This could result in the "mixing" of spectra with 
different amounts of absorption and normalizations, which 
cannot be well fit by the chosen model. 

Inspection of the fit residuals for the deep-dip spec- 
trum reveals clear structure around 6.4 keV, indicative of 
the presence of an Fe-K line. The addition of this line re- 
duces the x^ from 66.1 for 52 dof to 47.1 for 49 dof, which 
is significant at >99% confidence. The equivalent width of 
the line is 330 ± f 55 eV, the number of photons in the line 
is 4.4 X lO"** cm~^ s~^, the fine energy is 6.44 ±0.18 keV, 
and its width {a) is 310 ± 3^0 eV. The only parameter of 
the continuum model which changes significantly due to 
the addition of the line is the normalization of the COMPTT 
component, which decreases by 14%. There is no evidence 
for lines at 6.44 keV in the other dip spectra with la up- 
per limits to the number of line photons of 5.0x 10~^ cm~^ 
s-i and 7.1x10-^ 
dip intervals, respectively. 



5.2.2. Complex absorption model 

Since fitting the simplest model that provides statistically 
acceptable fits to the X 1658-298 dip spectra gives results 
which are inconsistent with the physical interpretation of 
the complex continuum model of Church & Baluciiiska- 



Church ( 1995 ), where the softer component is rapidly and 
strongly absorbed, we also fit this model to the same dip 
spectra. Since the blackbody and cutoff power-law model 
does not provide a good description of the persistent emis- 
sion spectrum (see Table ||), wc instead used the best- fit 
COMpTT and disk-blackbody model for this comparison. 
If the physical situation is similar, with the disk-blackbody 
originating from a compact region and the Comptonized 
component from an extended region, this should allow a 
meaningful comparison to be made with the complex con- 
tinuum model. 

In this case the model fiux is given by: 



,-ct7Vh 



e -" (/DBBe-"^« ± /cTT(/e-"^'« ± (1-/))), 
where /dbb and /ctt are the normalizations of the 
disk-blackbody and COMPTT components, / the cover- 
ing fraction of the COMPTT component, and a the photo- 
electric absorption cross-sections of Morrison & McCam- 



mon (1983). Note that the fundamental differences with 
the simple absorption model that is used in Sect. 5.2.1 are 



that the COMPTT component is partially absorbed and 
-^DBB and /ctt are fixed at the values measured during 
the persistent emission. 

The best-fit values of x^ are slightly worse than the fit 
results from the simple model, but show similar trends, 
e.g., the medium dip spectrum is not well modeled (prob- 
ably because of the same reason as discussed in Sect. 
|3.2.1 ) . The absorbing column for the disk-blackbody com- 



ponent smoothly increases from <9 x 10 atom cm" 
during the shallow dip to (3.73 ± 2.5) x lO^'^ atom cm~^ 
during the deep dip. The trends in the absorption col- 
umn and covering fraction of the COMPTT model are 
not as clear with the absorbing column increasing from 
(13.1 ± 2:2) ^ 10^^ atom cm"-^ during the shallow dip to 
>10^'^ atom cm^^ elsewhere. The covering fraction ranged 
from 0.31 ±0.02 during the medium dip to 0.91 ±0.01 dur- 
ing the deep dip. The results indicate that at all times the 
COMPTT component is obscured by a significantly higher 
column than the disk-blackbody. This confirms the same 
resul t obta ined using the simple absorption model (see 
Sect. ^.2.1 ) and is in contrast to the results previously ob- 
tained by applying the complex-continuum model to the 
dip spectra of other dipping sources. 
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Eclipse Spectrum 




ti^itr+V^ 



Energy (keV) 

Fig. 8. The LEGS and MEGS eclipse spectrum of X 1658- 
298 together with the best-fit absorbed power-law model 



5.3. Eclipse spectrum 

The spectral shape of the emission during eclipse was in- 
vestigated. Data from all four observed eclipses were com- 
bined to provide exposure times of 1.7 and 2.3 ks for the 
LEGS and MEGS, respectively. The HPGSPG and PDS 
spectra do not contain significant signals and were ex- 
cluded from this analysis. The resulting spectrum can be 
fit by an absorbed (Nh= (3.4 ± 1.6) x 10^^ atom cm"^) 
power-law with a = 2.0 ± H (Fig. |). The x^ is 35.3 for 
25 dof and the 2-10 keV luminosity is 4.3x 10^^ erg s~^ for 
a distance of 15 kpc. This is 2.7 ± 1.0% of the persistent 
value. Note the structured residuals in the lowest energy 
channels of the LEGS. This may suggest the presence of 
e.g., a dust scattered component, but the low count rate 
does not allow us to reliably investigate this further. We 
note that similar low-energy residual structure may also 
be present in the deep dip spectrum (cf. Figs, pi and 0). 
The 90% confidence upper limit to the number of photons 
in a 6.44 keV emission line is 7.6x10"^ cm~^ s""'^. 



6. X-ray bursts 

Nine X-ray bursts are visible in the light curve (see Fig. |l|). 
The first burst is the brightest with a peak flux measured 
during 2 s intervals of 1.5 times that of the other bursts. 
This burst was studied in detail to see if there was evi- 
dence for radius expansion which can be used to estimate 
the distance to the source. Light curves were produced 
in different energy bands. These show that some struc- 
ture is present in the burst, especially at higher energies 
(which is usually indicative of radius expansion). Time- 
sliced spectra (with 2 s time bins) of the burst where then 
obtained and the resulting spectra fit (after subtracting 
the persistent emission obtained in a 500 s interval be- 
fore the burst as background) with a blackbody model. 
However, the limited statistics do not allow us to clearly 
ascertain whether this burst shows radius expansion. The 



blackbody radius during the brightest 2 s of the burst is 
numerically larger than the other values, but is consis- 
tent (at the 1(7 level) with being constant. Using a flux 
estimate obtained from the two seconds when the burst 
is brightest (1.4 xlG~® erg cm~^ s""'^) and assuming that 
this corresponds to the theoretical Eddington limit for a 
1.4 Mq neutron star with a hydrogen atmosphere emit- 
ting isotropically a distance of 11 kpc is obtained. Using 
the empirically observed Eddington luminosity (3.0 xlO^^ 
erg s~^) of a sampl e of b ursters for which the distance is 
known (Lewin et al. 1993 ), a distance of 13 kpc is obtained. 
These estimates agree well with the distance of 15 kpc 



obtained by Gominsky (1981) who assumed that the 13 
burst sources observed by SAS-3 were uniformly dis- 
tributed around the galactic center. The distances to in- 
dividual burst sources, relative to the galactic center (as- 
sumed to be at 9 kpc), were obtained by comparing in- 
dividual properties with those of the sample. We have 
searched the literature for reports of other radius expan- 
sion bursts from X 1658-298; none were found. We con- 
clude that the luminosity of the flrst burst is likely to be 
at, or close to, the Eddington limit. Note that our esti- 
mates can be interpreted as upper limits to the distance, 
if the burst is not a radius expansion burst (and was thus 
emitting at below the Eddington luminosity). 

7. Discussion 

We present results of a 2000 August BeppoSAX observa- 
tion of X 1658-298. The persistent emission can best be 
modeled by the summation of disk-blackbody and Gomp- 
tonized components with low (1.3 x 10^^ atom cm~^) 
absorption. The best-fit to a 2-20 keV RXTE Propor- 
tional Gounter Array spectrum obtained 16 months ear- 
lier is to a cutoff power-law model with a = 2.1 ± 0.1, 
Nh= (5.0 ± 0.6) X 10^2 atom cm-^ and Ec = 8.6 ± 0.6 keV 
(WSB). When the same model is fit to the BeppoSAX NFI 
spectrum, the absorption is only 3.5 x 10^^ atom cm~^. A 
straightforward interpretation of this striking difference 
is that during the early stages of the outburst the X- 
ray emission is seen through a strongly thickened accre- 
tion disk, while by the time of the BeppoSAX observa- 
tion the disk had thinned and the central emission regions 
are directly viewed during non-dipping intervals. The ob- 
served 2-20 keV luminosity during the RXTE observation 
of 2.8x10^^ erg s~^ (WSB, converted from their quoted 
flux, which we assume to be the absorbed value) is indeed 
~40% higher than our 2-20 keV luminosity, consistent 
with a higher mass-accretion rate during the earlier phase 
of the outburst. 

The spectrum of the persistent emission from X 1658- 
298 can be described by disk-blackbody and Gomptonized 
components (the COMPTT model). However the x^ is not 
formally acceptable and low-level structure in the resid- 
uals is present which is probably due to the inadequate 
modeling of the source spectrum. The addition of Gaus- 
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sian line features improves the fit (without constraining 
the parameters of these features very well). This implies 
the presence of additional complexity in the spectrum of 
X 1658-298. The combination of a disk-blackbody or black- 
body model and a Comptonized component has been fit 
a large number of LMXRB spectra (see e.g., Guainazzi 
et al. |1998|; in 't Zand et al. |l999h| ; Barret et al. |2000| ; 



mal' component similar to that seen by Bonnet-Bidaud 
et al. (|200lD from EXO 0748-676 and arranged in a thin 



Sidoli et al. 2001), although a physical interpretation of 



the disk-blackbod y para meters is not straightforward (see 
e.g., Merloni et al. 2000). It is possible that the deviations 
from the smooth continuum in X 1658-298 are caused by 
intervening material, through which the emission regions 
are observed. The reasons why this intervening material is 
prominent in this source is probably due to the high incli- 
nation angle - the line of sight passes close to the surface 
of the accretion disk. Higher quality data from Chandra 
or XMM-Newton are needed to shed further light on the 
nature of these features. 

The spectra during dips can be reasonably well approx- 
imated with the persistent model, with only the COMPTT 
component undergoing strong additional absorption of 
^3 X 10^'^ atom cm~^. The normalizations of both spec- 
tral components decreases as the dipping gets deeper. The 
normalization of the disk-blackbody component shows the 
strongest dependence on total flux, decreasing to ~5% 
of the persistent value during the deep dip, compared 
to '^60% for the COMPTT component. This would sug- 
gest that the disk-blackbody component originates from a 
more compact region, while part of the Comptonized re- 
gion is still observed directly during dips. In this case it is 
likely that during the deepest dipping intervals the disk- 
blackbody is no longer viewed directly, the 5% of emisssion 
remaining could be X-rays scattered into our line of sight 
in an extended scattering region (the same region respon- 
sible for the COMPTT component perhaps) in a similar 
way to the 2.7% of emission still visible during eclipse. 
The compact nature of the disk-blackbody and the ex- 
tended nature of the Comptonizing region agrees with re- 
sults obtained by e. g.. Ch urch et al. ( 1997 ) and Church & 
Balucihska-Church ( 1995 ) with the "progressive covering" 
model. However, the strong additional absorption suffered 
by the COMPTT component is in contrast to the results of 
this model, where the (disk-)blackbody component rapidly 
undergoes strong additional absorption during dips. 

This additional COMPTT absorption is however consis- 
tent with the picture emerging from XMM-Newton obser- 
vations of the LMXRB EXO 0748-676, where the Comp- 
tonized component undergoes significant extra absorp- 



tion during dipping intervals (Bonnet-Bidaud et al. 2001), 
while the softer thermal component does not. However, 
during the EXO 0748-676 eclipses the thermal compo- 
nent remains, implying that it is the extended compo- 
nent in this system. It is possible to extend this picture to 
X 1658-298 to obtain a consistent picture of all the spectral 
changes observed during the dips. We first assume that the 
disk-blackbody component is actually an extended 'ther- 



layer above and below the accretion disk. In order that the 
normalization of the thermal component reduces to ~5% 
of its persistent emission value during deep dips, ~95% 
of the projected emission area should be obscured by the 
dense outer regions of the accretion disk. Assuming that 
the disk can be regarded as a dense solid body, then only a 
small fraction of the thermal flux will undergo absorption 
due to skin effects. This could explain why no additional 
absorption is detected from this component. Simulations 
of the structure of X-ray irradiated disks reveal that such 
sharp transitions in vertical structure may arise due to an 
ionization instability (Done & Nayakshin 2001). 

The Comptonized component does, however, suffer 
from additional absorption and from a ^40% reduction in 
normalization during deep dips. Assuming that this emis- 
sion originates from a spherical region centered on the 
neutron star, then during deep dips ^40% of the emission 
region should be obscured by the solid accretion disk as 
with the thermal component. The remaining ~60% is lo- 
cated away from the line of sight to the accretion disk. 
It is likely that the material in the line of sight to the 
Comptonized component is highly turbulent due to the 
interaction of the gas stream from the companion star 
with the outer accretion disk and the ionization instabil- 
ity described by Krolik et al. ( |l98l| ). This results in the 
formation of small relatively cool clouds embedded in a 
hotter medium (Frank et al. |l987|) . The passage of these 
clouds in front of the Comptonized region could result in 
the observed strong increase in absorption. 

While this model is totally ad-hoc, it is an attempt 
to explain the complexity observed during the dips. It 
does imply a particular geometry for the emitting regions, 
and in particular for the thermal component that must 
be located close to the accretion disk so that it can be 
almost totally obscured. We note that the two systems 
where this model may be applicable are the two eclipsing 
systems which are seen at inclinations angles of '^ 75-80°, 
compared to the pure dip sources which are viewed at in- 
cfinations angles of ^^60-75° (e.g., Frank et al. 1997). We 
speculate that the different viewing geometry may be one 
reason for the different spectral behavior during dips of 
X 1658-298 and EXO 0748-676 and the other dip sources. 

An interesting feature is the appearance of a strong 
Fe-K line during the X 1658-298 deep dip, whose energy 
is consistent with that of neutral or moderately ionized 
iron. Its origin (and absence during the other dip spectra) 
is rather difficult to understand, but maybe related to the 
amount of intervening material in the line of sight where 
substantial reprocessing of the continuum radiation can 
take place. 

The BeppoSAX ephemeris does not require a P term, 
in contrast to the ephemeris of WSB. This is probably 
due to the extremely uneven sampling (and possibly a 
small underestimate of the uncertainties) of the arrival 
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times in WSB, resulting in an artificial solution. We note Merloni, A., Fabian, A. C, & Ross, R. R. 2000, MNRAS, 313, 

that small (random?) deviations from a linear ephemeris 193 

are probably present, but this depends sensitively on the Morrison, D., & McCammon, D. 1983, ApJ, 270, 119 

assumed uncertainties in the arrival times. The observed Pamar, A. N., White, N. E., Giommi, R, & Gottwald, M. 

deviations may be similar in nature to those observed from t98b, Apj, Ms, 

EXO 0748-676 (see Hertz et al. ^. ^™ A^'lS' 3TO^'''''' °' ^' ^'^ ^'''"^'''' ^" '' ^^' ^^^^^ 
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